o 
o 

(N 



Modelling (001) surfaces of II- VI semiconductors 



m 






M. Ahr* M. Biehl T. Volkmann 

Institut fiir Theoretische Physik und Astrophysik 

Julius-Maximilians-Universitat Wiirzburg 

Am Hubland, 97074 Wiirzburg, Germany 

February 1, 2008 



cd 



-a 
c 
o 
o 



> 

o 

CO 

o 



cd 



c 
o 
o 



X 



Abstract 

First, we present a two-dimensional lattice gas model with anisotropic interactions which explains the exper- 
imentally observed transition from a dominant c(2 x 2) ordering of the CdTe(OOl) surface to a local (2 x 1) 
arrangement of the Cd atoms as an equilibrium phase transition. Its analysis by means of transfer-matrix 
and Monte Carlo techniques shows that the small energy difference of the competing reconstructions deter- 
mines to a large extent the nature of the different phases. Then, this lattice gas is extended to a model of a 
three-dimensional crystal which qualitatively reproduces many of the characteristic features of CdTe which 
have been observed during sublimation and atomic layer epitaxy. 



Within the last years, potential applications of 
electronic devices based on II- VI semiconductors have 
inspired basic research concerning for instance the 
surface reconstructions of these materials. For CdTe, 
a fairly complete qualitative phase diagram of (001) 
surfaces has been obtained ]jj]. An understanding of 
the interplay of these reconstructions with techniques 
like molecular beam epitaxy (MBE) or atomic layer 
epitaxy (ALE) is desirable both for technological ap- 
plications and for the theory of crystal growth. From 
the theoretical point of view, the case of CdTe is espe- 
cially interesting, since it exhibits a phase transition 
which involves competing reconstructions at typical 
MBE growth temperatures. While there have been 
a few models of MBE, especially for GaAs which 
take into account the effects of surface reconstruc- 
tions (e.g. Q), phase transitions between different 
reconstructions have not been considered in growth 
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models, so far. 

The zinc-blende lattice of CdTe consists of alter- 
nating layers of Cd and Te atoms the positions of 
which lie on regular square lattices. Under vacuum, 
the CdTe(OOl) surface is always Cd terminated. The 
surface is characterized by vacancy structures where 
the maximum Cd coverage pcd is 1/2. At low temper- 
ature, one finds a c(2 x 2) reconstruction, where the 
Cd atoms arrange in a checkerboard pattern. Above 
a critical temperature T c = 543XQ the surface is 
dominated by a local (2x1) ordering where the Cd 
atoms arrange in rows along the [110] direction which 
alternate with rows of vacancies jj], g, [|. Electron 
diffraction experiments have shown that there is a 
high degree of disorder in the (2 x 1) phase H. This 
suggests an interpretation of the reordering as an 
order-disorder phase transition. An external Cd flux 
stabilizes the c(2 x 2) reconstruction at T > T c . If a 
Te flux is applied, one finds surfaces terminated by 
Te dimers or trimers pL 0] . 



We first investigate whether the above mentioned 
transition can be explained within the framework of 
thermodynamic equilibrium. In a simplifying manner 
we treat the terminating Cd layer of the surface as 
a two-dimensional square lattice gas with sites either 
being occupied by a Cd atom or empty ||. Chem- 
ical bonds and the influence of the underlying crys- 
tal are accounted for by effective pairwise anisotropic 
interactions. An infinite repulsion excludes nearest 
neighbour (NN) pairs in [110]-direction (denoted as in- 
direction in the following). This can be justified from 
considerations like the electron counting rule [[[]. In 
[110](x)-direction an attractive interaction e x < fa- 
vors the occupation of NN sites. A competing attrac- 
tive interaction e,j < between diagonal neighbors 
(NNN) tends to stabilize the c(2 x 2) arrangement. 
The total energy of the system then reads 



H = e d n C d,d + £xncd,x - f-mcd, 



(1) 



where ncd, ncd,d and ncd,x denote the total number 
of Cd atoms, NNN and NN pairs in x-direction, re- 
spectively. The (effective) chemical potential \x con- 
trols the mean Cd coverage pcd — ( n x,y)- By choos- 
ing ed = — 1 the energy scale is fixed. In agreement 
with DFT calculations || , e x is chosen such that per- 
fect c(2 x 2) and (2 x 1) arrangements are nearly 
degenerate, the groundstate being a c(2 x 2) recon- 
struction. 

Besides the coverage pcd the system is character- 
ized by the two-point correlations C^ d = ncd,x/ncd 
and Cp d = ncd.d/ncd which determine the fraction 
of Cd atoms incorporated in locally (2x1)- or c(2 x 2)- 
ordered regions. 

Through transfer-matrix (TM) calculations and 
Monte Carlo (MC) simulations at constant pcd we 
obtain an estimate of the phase boundaries in the 
temperature-coverage plane (see fig. ||). At low 
temperatures (III), a c(2 x 2) ordered phase with 
Pcd ~ 1/2 coexists with a disordered phase of low 
coverage. At higher temperatures, the system be- 
comes homogeneously disordered (II) or ordered (I) 
depending on the coverage. For a range of coverages 
pcd the transition from (III) to (II) is accompanied 
by a sudden decrease of C^ d and a simultaneous in- 
crease of C% d , meaning that for these values of pcd 




Figure 1: Phase diagram of the 2D lattice gas with 
e x = —1.90. Solid lines and circles represent the 
phase boundaries, obtained from TM and MC. Right 
of the dashed line (squares) the c(2 x 2) structure is 
prevalent and vice versa. 



the (2 x l)-structure is locally prevailing in the dis- 
ordered regime (the area left of the dashed line (or 
squares) in fig. [j]). With decreasing \e x \ the tri- 
critical point, i.e. the point in the pcd~T -diagram 
where (I), (II) and (III) meet, shifts to smaller cover- 
age and higher temperature. Even more so does the 
line which separates c(2 x 2) from (2 x 1) prevalence. 
This might explain, why a transition from c(2 x 2) to 
(2 x 1) prevalence has not been found in ZnSe, so far. 
There, the energy difference between perfect c(2 x 2) 
and (2 x 1) reconstructions is indeed greater than for 
CdTe §. 

In order to investigate the evolution of the mor- 
phology of the surface in non-equilibrium processes 
like growth and sublimation, we extend the planar 
lattice gas to a model of a three-dimensional crys- 
tal. In a previous publication, this has been done for 
a simple cubic crystal which consists of alternating 
layers of the two atomic species M. In this paper, we 
present the next step towards a realistic modelling 
of CdTe by simulating the correct zinc-blende lattice 
structure. Besides the Cd-Te bond (binding energy 
e c ) we consider interactions between atoms on NN 
and NNN sites in one layer of the crystal. The basic 
idea here is to model the interactions of Cd atoms 
on the surface with the anisotropic interactions e x , 



ed and a hardcore repulsion in the y-direction while 
there is an isotropic attraction between Te atoms 
(tTe.b) and Cd atoms inside the bulk of the crystal 

Due to the topology of the zinc-blende lattice the 
condition of ergodicity can not be fulfilled within the 
limits of a solid-on-solid (SOS) model, where each 
atom must be bound to at least two atoms of the 
opposite species in the layer below. This can be 
achieved by considering the formation of Te trimers 
which have been observed experimentally Addition- 
ally, we consider Te atoms which have only one bond 
to a neighbouring Cd atom. In the following, Te 
atoms in these configurations will be denoted as Te* 
atoms. Te* atoms neutralize the hardcore repulsion 
between Cd atoms. Otherwise, the deposition of Te 
on a Cd terminated surface with a c(2 x 2) recon- 
struction would be impossible. 

As a simplification, we assume that point defects 
are not incorporated into the crystal. Then, the bulk 
can be described in a SOS manner, while the Te* 
atoms are considered seperately. In summary the 
Hamiltonian of the model is 

H = e c n c + ecd,bncd,b + tTe.b + ^h.ncd,h 

+ txiicd.x + edncd,d + £*n* (2) 

which is a function of the numbers of Cd-Te bonds 
(n c ), Te* atoms (n*), NN and NNN pairs of surface 
Cd atoms (ncd,x, ncd,d), NN pairs of Cd atoms in 
the bulk and Te atoms {ncd.b, nT e ,b) and surface 
Cd atoms next to incorporated Cd atoms (ncd,h)- 
Pairs of surface Cd atoms which are not bound to 
Te* atoms on neighbour sites in the ^-direction are 
forbidden. 

Besides the deposition of adatoms we consider des- 
orption, diffusion to NN and NNN sites and diffu- 
sion across steps. These are Arrhcnius-activated pro- 
cesses with rates vexp[—E/(kT)}. Dcsorption re- 
quires an activation energy E = AH. We have cho- 
sen Metropolis-like energy barriers E — max{B, B + 
AH} for diffusion processes. There is a barrier 
B = Bq for the diffusion of atoms with two bonds 
and a barrier B* for the diffusion of Te* atoms. 

Unfortunately, the available experimental data arc 
not sufficient for a systematic fit of the model pa- 
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Figure 2: pcd(+) and the correlations Cg d (x), 
C^ d (*) as functions of temperature in the stationary 
state of sublimation. The system size was 64 x 64. 



rameters and there are no values of diffusion barri- 
ers available from ab initio or semi-empirical calcu- 
lations. Therefore, the choice of the parameter set 
is partly based on estimates and physical reasoning. 
In the following, we present results that have been 
obtained with the parameter set e^ = — 0.086eU, 
e x = -0.168eV, e cd ,b = e Te .b = -0.069eU, e h = 0, 
e c = -0.43eV, e* = -0.26eU, B = 0.78eU, B* = 
O.UeV and v = 10 12 /s. 

Experiments have revealed a high density of steps 
on CdTe surfaces with a preferential orientation in 
the (100) direction. This leads to a dominant contri- 
bution of step flow sublimation. Therefore, we have 
simulated sublimation on vicinal surfaces with steps 
parallel to the [100] direction at a distance of 16 lat- 
tice constants. 

Our model reproduces the observation of Cd ter- 
minated surfaces in vacuum. A Te terminated sur- 
face decays in a temperature-dependent time inter- 
val! t(T). The corresponding desorption probabil- 
ity p(T) ~ 1/t(T) follows an Arrhenius law with 
an activation energy of 1.25 ± 0.02eV. After this 
onset, we obtain a stationary state where Cd and 
Te evaporate stoichiometrically. Then, the tempera- 
ture dependence of the sublimation rate follows an 
Arrhenius law with a higher activation energy of 
1.437 ±0.003eV. Experimentally, activation energies 




Figure 3: Correlations C^ d (dashed), C^ d (dotted) 
and thickness of the deposited layer (solid) in ALE 
as a function of time. The pulse time was Is, followed 
by 0.1s of dead time, the system size was 128 x 128. 



of 1.54eV |§ and 0.9eV || have been found for sto- 
ichiometric step flow sublimation and the decay of a 
Te terminated surface. 

Figure g shows pcd and the correlations C^ d , C^, d 
in the stationary state as functions of T. In the in- 
vestigated temperature interval, pea decreases only 
slightly from p cd = 0.49 at T = A50K to p cd = 0.43 
at T = 700K. At low temperature, the surface is 
dominated by a c(2 x 2) arrangement of the Cd atoms 
(C^ d > Cg d ). On the contrary, at high temperature 
the Cd atoms tend to arrange in the rows of the (2x1) 
reconstruction (C^ d < Cg d ). The crossover between 
both regimes is at T c = 57QK. This is the counter- 
part of the phase transition observed in the 2D lattice 
gas model. However, an investigation of C^ d , C^ d in 
the 2D lattice gas at the pcd measured in the 3D 
model shows, that the non-equilibrium conditions of 
sublimation increase the dominance of (2 x 1) over 
c(2 x 2) in the high temperature regime. 

ALE provides a scenario where our model can be 
investigated in a situation present in technical ap- 
plications. The basic idea is to obtain self-regulated 
layer-by-layer growth by alternate deposition of pure 
Cd and Te. Figure kl shows the evolution of the thick- 
ness of the deposited layer and the correlations C^ d , 
C% d at T = 500^ < T c . The initial state was a flat 



Te terminated surface. In the first phase, a Cd flux 
of 5ML/s is applied, in the second phase Te at the 
same rate. At the end of the first phase, the sur- 
face is Cd terminated with a c(2 x 2) reconstruction. 
In the Te phase, Te terminated islands form which, 
at the end of the cycle, cover « 50% of the surface. 
This is due to the fact that only half of the Cd atoms 
needed for a closed monolayer are present. Note, that 
the decrease of pcd at the beginning of the Te cycle 
leads to an increase of Cg d , which is consistent with 
the phase diagram (fig. 1) of the 2D lattice gas. In 
the following Cd phase, the islands are covered with 
Cd. Surprisingly, now the reconstruction is (2 x 1). 
This can be understood from the fact that the lattice 
gas interactions are present only between particles in 
the same layer. Thus, the island edges impose open 
boundary conditions to the lattice gas of Cd atoms 
on the island. In constrast to a c(2 x 2) reconstructed 
domain, a (2 x 1) terminated island can reduce the 
energy of its boundary by elongating in x-direction. 
Since the ground state energies of both structures are 
nearly degenerate, the formation of (2 x 1) may thus 
reduce the surface free energy. In the following Te 
phase, Cd and Te atoms diffuse into the gaps be- 
tween the islands, which leads to the formation of a 
closed, Te covered monolayer again. These observa- 
tions closely resemble ideas suggested in S to explain 
experimental observations. 

In summary, we have presented a model of 
CdTe(OOl) surfaces which reproduces several exper- 
imental observations in a semi-quantitative way. In 
the future, we will extend this model by consider- 
ing the dimerization of surface Te atoms and the fact 
that Te is adsorbed and desorbed as Te2- In cooper- 
ation with an experimental group, we will try a more 
quantitative determination of an appropriate param- 
eter set. 
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